Mon. Not. R. Astron. Soc. 000, ??-H5l (20091 Printed 7 August 2009 



(MN 1*1^ style file v2.2) 



On the Interstellar Medium and Star Formation 
Demographics of Galaxies in the Local Universe 

Matthew S. Bothwell^*, Robert C. Kennicutt, Jr^ and Janice C. Lee^'^ 

^ Institute of Astronomy, University of Cambridge, Cambridge, CBS OHA 
^ Carnegie Observatories, 813 Santa Barbara Street, Pasadena CA 91101 
^ Hubble Fellow 

Accepted 2009 July 30. Received 2009 July 17; in original form 2009 May 21 



ABSTRACT 

We present a demographic analysis of integrated star formation and gas properties for 
a sample of galaxies representative of the overall population at z 0. This research 
was undertaken in order to characterise the nature of star formation and interstellar 
medium (ISM) behaviour in the local universe, and test the extent to which global 
star formation rates can be seen as dependent on the interstellar gas content. Archival 
21 cm derived Hi data are compiled from the literature, and are combined with CO 
(J = 1 ^ 0) derived H2 masses to calculate and characterise the total gas content 
for a large sample of local galaxies. The distribution in stellar mass-normalised H i 
content is found to exhibit the noted characteristic transition at stellar masses of 
~ 3 X 10"'^*' Mq, turning off towards low values, but no such transition is observed in 
the equivalent distribution of molecular gas. Ha based star formation rates (SFRs) 
and specific star formation rates (SSFRs) are also compiled for a large (1110) sample 
of local galaxies. We confirm two transitions as found in previous work: a turnover 
towards low SFRs at high luminosities, indicative of the quenching of SF characteristic 
of the red sequence; and a broadening of the SF distribution in low-luminosity dwarf 
galaxies, again to extremely low SFRs of < 10^'^ Mq yr^^. However, a new finding is 
that while the upper luminosity transition is mirrored by the turn over in Hi content, 
suggesting that the low SFRs of the red sequence result from a lack of available gas 
supply, the transition towards a large spread of SFRs in the least luminous dwarf 
galaxies is not matched by a prominent increase in scatter in gas content. Possible 
mass-dependent quenching mechanisms are discussed, along with speculations that in 
low mass galaxies, the Ha luminosity may not faithfully trace the SFR. 
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1 INTRODUCTION 

The formation of stars is one of the primary drivers of galaxy 
evolution, and obtaining a full quantitative understanding of 
this process is an indispensable aspect of any investigation 
into galaxy behaviour. A profound symbiosis exists between 
the phenomenon of star formation and the form and content 
of the interstellar medium (ISM) of the host galaxy: a duel 
analysis, therefore, addressing both the star formation be- 
haviour and gas distribution in a sample of galaxies is ideally 
suited to characterising the properties of local galaxies, and 
analysing the results in the context of more global trends. 

Recent investigations into the evolutionary properties 
of galaxies have resulted in the striking finding that the dis- 
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tribution of galaxy colours exhibits a marked bimodality, 
separating 'red and dead' early type galaxies with little or 
no ongoing star formation, from a 'blue sequence' of star 
forming spirals , with the divide occu r ring at Mb ^ —19 
(Stratcva et al." '200l', 'B eU et al.l |2004 iBaldrv et all |2004 
ISchiminovich ct al. ,2007i) . This bivariate distribution can be 
examined in terms of many other galactic properties, such 
as the SFR; the colour-magnitude diagram can be inter- 
preted in terms of SF history and stellar mas s content (e.g. 
'Tromonti ot al." 2004, Br inchmann et al.ll20oi . iNoeske et all 
(2007 , .Salim et al. 2007 ). 

The divide between the red and blue sequences emerges 
also in the star formation rate per unit stellar mass (Spe- 
cific Star Format i on Ra te, S SFR), as demonstr a ted b y 
iBrinchmann et al.l and ISchiminovich et all l|2007h : 

the distribution in the SSFR vs. M, plane is bimodal, 
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with a red sequence of low SFR massive galaxies (SSFR 
~ 10~^^ yr~^), and a less massive blue sequence of galax- 
ies with significatly high er SSFRs (10"^ - 10"^° yr"^). 
ISchiminovich et alj |2003) use these areas to classify galax- 
ies as either on or off a 'SF sequence', along which many SF 
parameters (such as the SFR surface density, Esfr) are close 
to constant. In this picture, the blue 'main sequence' galaxies 
evolve secularly towards higher M, while gradually lowering 
their SSFRs, while the red sequence turn over is caused by 
the quenching of the SF. Conversely, galaxies can be driven 
off the main sequence in the other direction (towards higher 
SSFRs) by an event - such as an interaction - that triggers an 
increased SFR, or a sudden starburst (e. g. Ifernauistlll989l : 
iBarnes fc Hernauistlll99ll : iHopkins et allbood )^ 

The rate at which a galaxy is forming stars is, however, 
in some sense a 'high level process', being both complex and 
dependent on the more fundamental underlying properties 
of the galaxy - primarily the content of the ISM. The galac- 
tic 'gas fraction' (as opposed to the absolute amount of gas, 
which is strongly dependent on the mass of the host galaxy) 
provides information on the evolutionary history of the sys- 
tem in question, and also determines the potential for future 
star formation. As such, an understanding of how the gas 
fraction varies from galaxy to galaxy may help shed light on 
the origin of the bifurcation in the colour-magnitude plane; 
tentative explorations of the causes underlyin g this phenom- 
ena h ave highlighted the role of gas - iDekel fc BirnboimI 
l|2006l ). for example, note that the bifurcation mass coin- 
cides with the transition between the 'cold accretion' and 
'cooling flow' regimes of model galaxy formation. However, 
the role of interstellar gas in the development of the range 
of galaxies we see today has not been fully explored. A key 
question to address is whether the bimodal division into two 
galactic types can ultimately be seen as a consequence of the 
underlying gas behaviour. 

Investigations into the global properties of galaxies have 
principally focused on the upper end of the luminosity func- 
tion, due to the ease at which brighter systems can be stud- 
ied. However, investigations into star forming behaviour in 
low mass dwarf gala xies have resulte d in some unexpected 
findings. A study by iLee et al.l (|2007l ') investigated the dis- 
tribution of the SSFR (measured via the Ha Equivalent 
Width) with absolute magnitude for a small volume limited 
sample of local galaxies (obtained from the llMpc Ha UV 
Galaxy Survey, [IIHUGS]), with Ha derived SFRs for the 
population of local low luminosity dwarf galaxies, down to 
M b ~ —10. In addi tion to the well studied transition (found 
bv lLee et al.) (|2007l ) to be at Mb = -19), a second transition 
was found, with low mass dwarf galaxies (Mb > -15) having 
a very large (> 3 dex) spread in their specific star formation 
rate (SSFR). Between the two transition regions, the SSFR 
was found to inhabit a very restricted range of values, imply- 
ing a high degree of regulation in the star forming behaviour 
of the blue sequence - this comprises the 'main sequence' of 
star forming galaxies discussed above. These findings sug- 
gest that there may be three distinct modes of star forming 
behaviour operating across the mass spectrum, rather than 
two as previously thought. A lack of cor relat ion with envi- 
ronmental factors (|Hunter fc Elmegreenll2()04h suggests that 
internal processes are responsible for this third mode, but 
the nature of these process is, as yet, unclear. 

In this paper, the star formation properties and inter- 



stellar medium (ISM) content of a sample of local galaxies 
are characterised in order to delineate the various popula- 
tions that exist along the low-redshift Hubble sequence. As 
such, we seek to determine the extent to which the previ- 
ously described bimodality in the SSFR is solely a result of 
the s tandard scaling of star formation rate with gas content 
(e.g. ISchmidt|[l959l : iKennicutjl 1998bl ) . 

Our primary sample, the IIHUGS dataset, is highly 
complete at the lower end of the mass function, but suffers 
from a paucity of bright massive objects: while over 80% of 
the sample has a SFR lower than the LMC, just ~6% have 
SFRs greater than th e Milky Way (taking SFRm w=2-6 Mq 
yr~^, as reported by iBoissier fc Prantzoj|l999l ). Therefore, 
we aim to augment the IIHUGS dataset with more volumi- 
nous surveys, to see whether the inclusion of more data at 
the brighter end of the luminosity function would corrobo- 
rate previous findings. Our overall compiled sample, consist- 
ing of nearly a thousand local galaxies, provides a dataset 
with not only excellent completeness at low (~SMC) masses, 
but also coverage over a large enough volume to statistically 
sample galaxies with masses comparable to or greater than 
the Milky Way. This sample allows star formation rates and 
gas contents to be compared in a uniform manner over a 
mass spectrum spanning five orders of magnitude. This pa- 
per will also re-evaluate the possibility of a 'third mode' 
of star formation existing in low mass dwarfs, again in the 
context of the available gas supply. 

In §2, we describe the samples and data drawn upon for 
this work, and briefly describe the calculations we have used 
to derive various parameters. In §3, a demographic analysis 
of the gas behaviour, star forming behaviour, and the link 
between the two are presented in turn for both a highly 
complete local volume sample, and a more populous deeper 
sample. These results are discussed in more detail in §4, and 
we present conclusions in §5. 



2 METHODOLOGY 
2.1 Sample Description 

There were two primary drivers behind the selection of data. 
Firstly, a large sample was required, in order to represent as 
fully as possible the trends in gas properties and star form- 
ing behaviour in the local universe. In order for this aim to 
be realised, the sample must be relatively homogeneous, in 
order to investigate the range of 'normal' behaviours; many 
galaxies are observed precisely because of their unusual be- 
haviour, with the result that an uncritical compilation of Ha 
and H I data from the literature would likely result in biases 
towards more extreme objects. Secondly, in order to char- 
acterise the full range of behaviours across the entire mass 
spectrum, good completeness at the lower end of the lumi- 
nosity function was required. Low mass dwarf galaxies are 
the numerically dominant galaxy type, but Malmquist-type 
biases result in their exclusion from many observational in- 
vestigations. As a result, volume limited local surveys (which 
will be dwarf-dominated) have been utilised. The resulting 
combined sample has the advantage of being highly com- 
plete at low masses, and large enough to be representative 
at high masses. In order to ensure that the final compilation 
is as bias-free as possible, the analysis has been repeated for 
the local volume sample (see §3.2), with no change in result. 
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2.2 Hq 

Information on the star formation properties of galaxies 
in the Local Vo l ume w as obtained from the Hq survey of 
iKennicutt et al.l (|2008l ). This survey is an essentially vol- 
ume limited sample of 436 local galaxies (D ^ 11 Mpc). 
Galaxies were selected with the added proviso that galaxies 
in the so called 'zone of avoidance' (e.g. lying in the galactic 
plane, |b| > 20°) were excluded. This survey includes the 436 
known galaxies within 11 Mpc, and is essentially complete 
for Mb ~ -15.5, and M(Hl) ~ 10** Mq (Lee et al. 20(3), 
providing excellent statistics on the star formation rates of 
the low luminosity dwarf population. 

This sample has good statistical sampling of the lower 
end of the luminosity function, but due to the rarity of bright 
(Mb < -19) galaxies in the local volume surveys covering 
larger volumes were required to achieve improved sta t istica l 
coverage at the upper end. The survey of lJames et al.l (|2004l ) 
was chosen for this purpose. This is a Hq based survey of 334 
galaxies, taken from the Uppsala General Catalogue (UGC), 
which has a magnitude limit of B = 15.5. The galaxies ex- 
amined in the James et al. Hq survey were chosen to have 
recession velocities of < 3000 km s~^ (D ^ 43 Mpc; assum- 
ing Ho = 70 km s~^Mpc~'^). James et al. also stipulated 
that the galaxies have a morphological T-type ^ (SO/a or 
later) . 

In addition to the sample of lJames et al.l l|2004l ). many 
galaxies were included from the large database of Hq galax- 
ies (Kennicutt et al. 2010, in prep). More than 3000 galaxies 
have been observed in Hq l|Kennicutt et ah I l2008l ). provid- 
ing a rich and varied dataset to sample. However, the data 
are extremely heterogeneous, with many of the galaxies be- 
ing investigated precisely because of their unusual, or indeed 
exceptional, behaviour. In order to provide a characteristic 
representation of local universe star formation, unbiased to- 
wards extreme objects, surveys and compilations specifically 
ta rgeting cluster mem bers (such as the GOLDMine sample 
of lGavazzi et al.ll2003l ) and starburst galaxies were excluded. 
The complete Hq sample consists of 1110 galaxies. 

2.3 HI 

Detection of the neutral atomic component of the ISM in low 
redshift galaxies is comparatively easy, d ue to the high sen- 
sitivi ty available at wavelengths of 21cm (jRoberts fc Havned 
1 19941 ) , and as such many large H I surveys exist in the liter- 
ature. For this study, information on the neutral hydrogen 
content of galaxies was obtain ed from the la r ge ho moge- 
neous compilation presented bv lSpringob et al.l(|2005h . This 
is a large (~ 9000) galaxy compilation of optically selected 
galaxies in the low redshift universe, collected over a pe- 
riod of 20-1- years, selected to have heliocentric velocities of 
-200 < V0 < 28, 000 kms"\ A total of 8850 galaxies were 
detected unequivocally in Hi with a further 156 and 494 
being marginally detected and not detected respectively. 

In order to remove the strong effect of galactic stel- 
lar mass on the Hi content, and to provide a parameter as 
closely analogous to SSFR as possible, the Hi content is 
presented in the form of a Hi mass to stellar mass ratio, 
7?(Hl), defined as Mhi/M*. This constrains the analysis to 
those galaxies with available photometry in order to measure 
the stellar mass (see §2.5 below). A total of 7486 'detected' 




Figure 1. Schematic Venn diagram illustrating the 3 data types, 
respective sample sizes, and number of galaxies with combinations 
of data. Note that the numbers in the coloured segments represent 
the total number of galaxies recorded with the respective data 
type, as opposed to the standard Venn formalism whereby the 
numbers represent galaxies recorded exclusively in the respective 
data type. 

galaxies from the Springob survey possess archival B-band 
photometry - this makes up the bulk of the Hi sample. 

In addition, to provide Hi information for loca l dwarf 
galaxies, the H I data from the iKarachentsev et al.l (|2004l ) 
compilation will be included in the analysis. This is an all 
sky catalogue of 451 galaxies, selected to have a distance 
^ 10 Mpc (= Vrad ^ 550 km s~^), which contains both 
optical photometry and 21cm Hi data, making it ideal for 
an analysis of ISM content. The limiting magnitude of the 
survey was Bt = 17.5, with a resulting absolute magnitude 
limit of Mb ~ -7 at 1 Mpc, and Mb ~ -12 at 10 Mpc. 

Additional 21 cm Hi data to complement the Ha 
dataset was obtained from the sample of Springob et al., 
the NASA Extragalactic Database (NED), and the on- 
line reference catalo gue of astrophysical data HYPERLEDA 
(|Paturel et al.ll2003l ). 

The total number of galaxies included in the H I analysis 
is therefore 7819. 

2.4 Ha 

The study of molecular hydrogen is important to the under- 
standing of star formation: to the best of our knowledge, it 
is only inside the giant molecular cloud (GMC) complexes 
that stars can form. H2 is much more challenging to de- 
tect than Hi; few galaxies hav e been mapp e d in CO (the 
FCRAO extragalactic survey of I Young et ahl l|l995l ), for ex- 
ample, contains resolved CO data for 300 galaxies), and 
fewer still have the unresolved CO measurements required 
to measure the global CO flux, necessary to calculate a total 
molecular hydrogen mass. 229 galaxies in the FCRAO sur- 
vey have global integrated CO fluxes recorded. In addition, 
archival CO data for 245 galaxies has been compiled and pre- 
sented bv lObreschkow fc Rawlingd (|2009l ) , which were added 
to the sample. H2 is very difficult to detect in dwarf galax- 
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ies; they have httle to no CO emission, and uncertainties in 
the CO-H2 conversion factor at low metallicities mean that 
it is extremely challenging to get accurate H2 masses for the 
lowest luminosity galaxies in our sample. A recent survey by 
iLerov et al.1 (|2005l ) attempts to address this issue, and from 
there we take 10 sources that overlap with our sample. A 
total of 392 unique galaxies with single dish CO fluxes were 
used, 43 of which have only upper limits on their CO fluxes. 

Figure [T] shows a diagrammatical representation of the 
sample sizes for different data types, and the number of 
sources with various combinations of available data. 



2.5 Derived Parameters 
HI mass: 

Hi, detected via the hyperfine splitting spin transition of 
atomic hydrogen emitting at 21 cm [v — 1420 MHz), is the 
most powerful available probe of the ISM. The forbidden 
nature of the transition results in most accumulations of H I 
being optically thin, with a direct proportionality between 
the 21 cm luminosity and the number of emitting atoms. 
The neutral hydrogen mass is therefore given by the simple 
formula 

Mhi = 2.36 X lO'' J Sdv Mq 

Where J Sdv is the integrated 21cm line flux in Jy km s~^, 
and the distance D is in Mpc. 

H2 mass: 

The molecule H2 lacks a permanent electric dipole mo- 
ment, and as such is very difficult to detect directly. In- 
stead, CO emission (e.g. the J = 1 ^ transition at 
A = 2.6mm) is measured, and converted to a H2 column 
density (e.g. lYoung et al.lll995l ) The 'standard' calibration 
used (the oft discussed CO-to-H2 conversion factor, given by 
[Kenn cv 1997) is 

Nh^/Ico = Xco X 10^" H2 cm-2 [K(Tr) kms-i]-i 

Where Xco is generally taken to be 2 - 3; we take the X- 
factor to be a constant value of 2.8. Kenney and Young 
(1989) show that the molecular hydrogen mass can therefore 
given by 

Mh2 = 1.1 X 10* ■ Sco Mq 

Where Sco is the global J = (1 ^ 0) flux in Jy km s"'^, and 
the distance D is in Mpc. This assumption of a constant CO- 
to-H2 conversion factor (as oppos ed to, say, a luminosity 
depe ndent factor as presented by lObreschkow fc Rawlingg 
I2OO9I ) will have the tendency to underestimate the mass of 
molecular hydrogen in the low mass dwarf galaxies. However, 
the quantities involved are negligible fractions of the total 
mass of the ISM; for this analysis, therefore, the constant 
factor has been used. 

Absolute Magnitudes: 

Absolute magnitudes (Mb) were derived from the total B- 
band magnitudes corrected for intern al extinction and incli- 
nation effects (Bt) given in the RC3 of lde Vaucouleurs et al.l 
l|l99l[l , using the standard conversion formula 



M = m - 5 log Dl + 5 

Where corrected B-band magnitudes were unavailable, the 
formulae given by de Vaucouleurs et al. were used to convert 
total B magnitudes into this system: 

Bt = Bt - Ag - Ai - Kb 

Where Bt is the total B-band magnitude, Ag is the 
galactic extinction (calculated here using the IR maps of 
ISchlegel et alll 1 998*1. Kb is the K correction for redshift, and 
Ai is a linear function of log R25 (which here is the isophotal 
axis ratio, as opposed to an isophotally-defined radius), 

Ai = a(T)logR25 

with 

a(T) =0 T < 

q(T) = 1.5-0.03(T-5)^ T^O 

The K-correction has been taken to be negligible, as the 
most distant galaxy is UGC 12840 at 93 Mpc {z = 0.023 for 
a concordance cosmology of flm = 0.3, 57a = 0.7, h — 0.73), 
and the vast majority are closer than 40 Mpc {z < 0.01). 

Stellar Mass: 

iBell fc de Jo"^ l|200ll ) note that the stellar mass to light ra- 
tio shows trends with various galaxy parameters (including 
surface brightness, gas fraction, K-band magnitude, colour, 
and galaxy mass), and - most critically for this study - has 
a dependency on the recent star formation history of the 
galaxy in question. Using any given photometric luminos- 
ity (even the K-band) as a proxy for the stellar mass will 
therefore introduce errors, as there is no passband where 
the mass to light (M/L) ratio is constant. Bell & de Jong's 
model (M/L) ratios correlate strongly with galactic colour; 
they therefore provide the algorithm 

log(M/L) — a\ + bx ■ colour 

where for B - V colours and B-band photometry, axCB) = 
-1.110 and &a(B) = 2.018, for a scaled Salpeter IMF and a 
mass dep endent formation epoch model (as described by, for 
example, ICole et al.|[200ol ) . 

A potential problem with this method is the lack of 
complete photometric datasets from which the colours can 
be derived: archival B - V photometry (corrected for redden- 
ing) is only available for ~ 50% of the sample (498 galaxies of 
the total 1049). However, it is possible to ac curately estimate 
a colour based on morp hological type (e.g. lHolmber^ll958l : 
[Roberts fc HavIle^ll994^ . where the morphological types are 
coded in the RC3 system as reported in the Hyperleda 
datab ase. The colours are taken from [Roberts fc Havnea 
(|l994l ). and are shown in Table [T] Figure [2] shows the stellar 
mass derived from the real and assumed colours for the 498 
galaxies with B - V photometry. It was found that using the 
colours derived solely from morphological type data caused 
the stellar mass to be significantly underestimated for the 
lowest mass dwarfs. To correct for this effect, galaxies with 
masses calculated to be < 10^'^ Mq were assigned a (B- 
V) colour of 0.3 - this resolved the issue. The two methods 
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lo*" 10' lo** 10' 10'° 10" 10'^ 

stellar mass: {B-V)^^^^,^! [solar masses] 

Figure 2. Plot of the stellar mass, derived using the method of 
Bell and de Jong (2001). The y-axis displays the stellar mass de- 
rived using real B-V photometric colours, and the x-axis displays 
the stellar mass derived using the assumed B-V colour given in 
Table n 

Table 1. Morphology-Colour Relationship, with resultant B- 
band mass-light ratio (M/L): colours from Roberts & Haynes 
(1994), M/L ratios calculated using the algorithm of Bell & de 
Jong (2001). *Note that galaxies with stellar masses < lO'"'^ Mq 
were assigned a (B-V) colour of 0.3 - see text for details. 



T-type 


< B- V > 


M/L 


SO 


0.9 


5.08 


SOa 


0.8 


3.19 


Sa 


0.7 


2.00 


Sab 


0.7 


2.00 


Sb 


0.6 


1.26 


Sbc 


0.55 


1.00 


Sc 


0.5 


0.79 


Scd 


0.5 


0.79 


Sd 


0.5 


0.79 


Sm 


0.4 


0.49 


Im 


0.4 


0.49 


dirr* 


0.4 


0.49 



for mass estimation agree well, with reasonably small scat- 
ter (see i)4.1l for discussion) and no significant trend in the 
residuals: this suggests that the colour assumption method 
is a acceptable one for the estimation of stellar masses; this 
assumptive technique will be used to estimate stellar masses 
for the remaining 40% of the sample for which no colour is 
available. It is important to note that restricting the analy- 
sis in this work to only those galaxies with archival colours 
(and, therefore, reliable stellar mass estimates) does not al- 
ter the scatter or trends in any of the derived relationships, 
and does not alter any of the conclusions drawn. 

SFR: 

Photoionisation of interstellar gas by massive (> 10 Mq) 



stars results in Balmer line emission, with a flux propor- 
tional to the number of incident ionising UV photons. By 
measuring the global flux in a Balmer emission line (no- 
tably Ha), it is possible to calculate the ionising UV flux, 
and therefore the number of O and B stars in a galaxy. O 
and B stars being very young, it is possible (assuming some 
IMF) to extrapolate to the total instantaneous star forma- 
tion rate. iKenniciittI l|l998ah gives the conversion 

SFR (M0 year)-' = 7.9 X 10""^ L(Ha) (ergss"') 

Which is based on a standard Salpeter IMF (a = 2.35) from 
0.1 - 100 Mq. Most narrowband surveys for Ha use filters 
which are wide enough to detect the adjacent [NII]6548,6583 
lines, and simply publish the combined Hq:-|-[NII] flux. In 
such C&S6 S, cl correction for [Nil] contamination must be ap- 
plied (see lKennicutt et al.|[200a ): 

log([NII]/HQ) = (-0.173 ± 0.007)Mb - (3.903 ± -0.137) 

for Mb > -21, and 

log([NII]/Ha) = 0.54 

for Mb ^ -21. 

2.6 Nebular Attenuation Corrections 

The internal Ha extinction is crucial to constrain in order 
to accurately determine star formation rates. Internal ex- 
tinction due to dust is the largest source of systematic error 
when measuring SFRs (|Kennicuttll998ah . and the large vari- 
ations between galaxies can also introduce signiflcant ran- 
dom errors. 

The simplest (and crudest) method of correcting for ex- 
tinction is to adopt a constant e xtinct ion factor (1.1 mag- 
nitudes was used by iKennicuttI (|l983l ). for example); this 
method is lacking in rigour, however, and is insufficient for 
a thorough analysis, as the r eal amount of ex: tinction can 
vary widely between galaxies l|Buat fc Xull 19961 ) . A method 
which takes into account the general trend that the attenu- 
ation increases with increasing galaxy luminosity (i.e. dwarf 
galaxies are more trans parent on a verage than more massive 
spiral discs) is given in ILee et al.l (|200^) as: 

A(Ha) = 0.14 (Mb > -15) 

A(Ha) = 1.971 + 0.323Mb + 0.0134M| (Mb < -15) 

This is calculated from the Balmer decrement, using 
the sp ectroscopic measurements of lMoustakas fc KennicuttI 
(|2006l ). A more rigourous method is to calculate the Ha 
extinction directly from the total infrared (TIR) luminos- 
ity of a galaxy. This method is based on an 'energy bal- 
ance' assumption that the extincted Hq flux will be re- 
radiated in the IR (se e iKennicutt et al. 1 12009| for details). 
IKennicutt et al.l (|2009l ) give the prescription 



A(Ha) =2.5 log 14-0.0025 




for calculating A(Ha), the extinction in Ha. Due to the 
self-similar shape of the infrared spectral energy distribution 
(SED) at different luminosities, the total IR luminosity can 
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be accurately estimated (to better than 1% for most galax- 
ies) using a superposition of 2 5A(m, 60^m, and 100/ im IRAS 
fluxes, using the algorithm of lOale fc Heloul (j2002l i: 



L(TIR) = Cii^L„(25Mm) + C2i^L^(60/im) + Cs,vh^(im^Jiin) 

where Ci = 2.403, C2 = -0.2454, and Cs = 1.6381 at 
z = Q. Estimating the extinction in this way does require 
that the galaxy in question has been observed by IRAS, 
which is only true of a small percentage of the total sample. 
If no IRAS fluxes are available, the extinction will be taken 
to be the quadratic function of Mb described above. 



3 ANALYSIS 

In order to investigate the properties of our sample of galax- 
ies, a demographic analysis has been carried out on two char- 
acteristic parameters, i) _R(gas) and ii) the star formation 
rate (both of which are closely linked via the Schmidt law 
formalism), and their trends with absolute magnitude. 



3.1 Gas Demographics 

In this analysis of the gas demographics, the atomic and 
molecular pha ses of the ISM a r e pres ented separately. First, 
the sample of ISpringob et aL* (2005") is combined with the 
Local Volume Hi data of Karachcntscv ct al. (2004), to pro- 
vide a large compilation of _R(Hl) {= Mhi/M*). Secondly, 
the available H2 data are presented in the same way. Finally, 
the H I data are combined with H2 data to present total gas 
ratios, 7?(total). 



3.1.1 Atomic Hydrogen 

In order to remove the effect of galaxy mass on the H I con- 
tent, we normalise to the stellar mass, and present _R(Hl), 
defined as Mhi/M«. Figure |3] plots -R(Hi) against absolute 
B-band magnitude for the entire Hi sample as described 
above. 

Taking the Hubble sequence to be a rough function of 
luminosity, it is clear t hat Figure [3] reproduces the (by now) 
well known result (e.g. [Roberts fc Haviie3ll994h that the H I 
ratio increases along the Hubble sequence, with the smaller 
(late type) galaxies typically having7?(Hl) > 0.5 (and most 
having i?(Hl) > 1), and the mean value of i?(Hl) decreasing 
as luminosity increases (and the population shifts to being 
made up of earlier types). The dense centre of the plot (i.e. 
the peak of the sample's number distribution) is dominated 
by intermediate type (Sa, Sb, Sc) discs, which have typical 
7?(Hl) of ~ 0.2. 

The distribution is smooth for Mb > —18, with the 
mean gas mass ratio decreasing monotonically with increas- 
ing luminosity. At Mb < —18, however, the distribution is 
joined by outliers existing below the main locus, as early 
types appear with extremally low iJ(Hl) of < 0.05. The en- 
tire distribution (in log-T? - log-Ls space) can be fit by a 
single power law of index —0.58 ± 0.007 (see i ]3.2.2l for a 
fuller discussion of the slopes of the distributions). 
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Figure 3. -R(Hl) plotted against absolute B-band magnitude for 
the 7819 galaxies from the H I sample with concomitant B-band 
photometry and morphological type data. Contours have been 
used where the density of points is high; where points become too 
rarified for contours to faithfully trace the distribution, individual 
points have been plotted. The Least-squares best fitting power law 
fitted to the distribution has an index of (—0.58 ± 0.007). 



3.1.2 Total Gas Content 

The cool, neutral phase of the ISM, composed of atomic 
hydrogen, acts both as a reservoir of fuel for future star 
formation, and as a tracer of star formation (as new massive 
stars ph otodissociate the surr ounding molecular ISM into 
Hi - e.g. lAndersson et al.lll99ll ). Stars form inside molecular 
clouds, however, which are comprised primarily of molecular 
hydrogen - H2 - and helium, and thus an understanding of 
both the neutral and molecular components of the ISM is 
integral to a full characterisation of the relationship between 
gas content and SFR. To account for the mass of helium, a 
standard correction is applied with most authors taking the 
helium mass to be a fixed fraction of the total (Hi -j- H2) 
hydrogen mass, i.e. 



Mgas = 



(Mm -fMH,) 
0.74 



Figure |4] shows the relationship between i?(gas) and 
B-band magnitude for the 392 galaxies with both Hi and 
CO-derived H2 data. The upper, middle, and lower panels 
show the Hi mass ratio, H2 mass ratio, and total mass ra- 
tio respectively, all plotted against Mb . Galaxies with upper 
limits on their CO fluxes (10% of our sample, as described 
in §2.4) have been marked with arrows in panels (b) and (c). 

Plot (a) demonstrates very similar behaviour to the 
larger H I sample shown in Figure |3] (This new set is highly 
incomplete at the less massive end, with only 5 galaxies with 
Mb ^1 — 15 having recorded H2 masses). Late morphological 
types have, on the whole, larger gas mass ratios than their 
early-type kin, making up the purple and blue points on the 
upper side of the distribution. As with the larger sample, at 
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Figure 4. Upper panel: iJ(Hl) against Mb for all 392 galax- 
ies with recorded Hi and H2 data. Middle panel: -R(H2) plotted 
against Mb for the same sample. Lower panel: _R(total), (defined 
in text) plotted against Mb for the same sample. A dashed hori- 
zontal line has been drawn at a gas mass ratio of 50%, for ease of 
reference. Galaxies with only upper limits on their CO mass have 
been marked with arrows in panels b and c. 



the transition luminosity Mb = — 19 the distribution does 
show a population of gas-poor outhers, with some galaxies 
brighter than this having a marked decrease in their neutral 
gas mass ratios, down to extremely low values of < 0.5% . 

Plot (b) shows -R(H2) for the same sample. The dis- 
tribution in H2 is markedly different to that in Hi: the plot 
shows no significant trend in -R(H2) with Mb, supporting the 
well known idea that the amount of molecular hydrogen in 
the ISM of low mass galaxies is much lower than their more 
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Figure 5. Star Formation Rate plotted against absolute B-band 
magnitude for the complete Ha sample. The solid diagonal line 
represents a value of P* = 0, while the upper and lower dashed 
lines show P* = 1 and P* = — 1 respectively. The gray vertical 
lines at Mb = -15 and Mb = -19 indicate the boundaries sep- 
arating the three SSFR behavioural regions identified by Lee et 
al. (2007). Galaxies with only upper limits on their H« fiuxes are 
marked with downward arrows. 



massive counterparts (or equally that it becomes harder to 
detect due to variations in the H2-CO relationship). There 
also is an apparent limit, at a molecular gas mass ratio of 
^ 50% (log i?(H2) = 0), which exists independent of the 
mass of the host galaxy. 

Plot (c) shows _R (total) for the same 392 galaxies. In 
less massive galaxies, plot (c) is broadly similar to plot (a), 
albeit with (very slightly) larger values of _R(gas). These late 
type galaxies are dominated by neutral hydrogen, and do 
not have a significant molecular component to their ISM. At 
the upper end, however, the molecular contribution becomes 
significant. The outliers, existing towards lower gas mass 
ratios seen in the H I plots, are much less apparent, in favour 
of a smoother, more monotonic decrease. 



3.2 Comparison of Gas and SFR Demographics 

In this analysis of the SFR demographics, the sample of Ha 
detected galaxies has been divided into two parts. Firstly, 
the entire Ha sample - comprising 1110 galaxies - has been 
analysed, in order to investigate trends as completely as pos- 
sible. Secondly, in order to directly compare the SFR prop- 
erties with the gas content, the sample has been restricted 
to only those galaxies with both Ha and H I data. This com- 
posite sample comprises 1049 galaxies (see Figure [T] for clar- 
ification) . 

Figure [5] shows the SFR plotted against absolute B- 
band magnitude for all galaxies with available Ha data. The 
vertical dashed lines represent the two transitional bound- 
aries previously identified. The diagonal lines represent dif- 
ferent values of the parameter P*, defined bv lKaisin et al.l 
|2003) as 



P* = lot 



/ SFR - To 
V Lb 
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Where Tp is the Hubble t ime, taken here to be 13.7 x 10^ 
years l|Spergel et alj 12007). This par ameter is analogous to 
the Scalo 'b' parameter jScalolll986l ). in that a galaxy with 
P* = is forming stars at such a rate that if it were to form 
stars at the current rate over one Hubble time, the total 
produced luminosity would equal the current galactic lumi- 
nosity. By extension, values of P* = — 1 and P* = 1 imply 
that it would take times of 10 To and 0.1 To respectively to 
fulfil the same condition. 

The central solid line in Figure [5] represents a value of 
P* = 0, whereas the upper and lower (dot-dashed) lines 
represent P* = 1 and P* = — 1 respectively. With the vast 
majority of galaxies lying between these upper and lower 
limits, it is possible to define a 'main sequence' of star form- 
ing galaxies this way, with (1 > P* > —1). As expected, 
the mean value of P* decreases with increasing luminosity. 
Conspicuously, however, there are many points existing as 
extreme outliers, with SFRs far below this 'main sequence'. 
These occur almost exclusively in the upper (Mb < —19) 
and lower (Mb > —15) luminosity categories, with very few 
galaxies in the intermediate bin having aberrant SFRs. The 
three previously identified modes of SF behaviour are evi- 
dent in Figure [5] A high mass population containing some 
galaxies with quenched SF, which belong to the 'red se- 
quence'; a secularly evolving intermediate population, with 
highly regulated SFRs; and a low mass, low luminosity pop- 
ulation, again containing members with anomalously low 
SFRs. Within the 'main sequence' (i.e. excluding those aber- 
rant galaxies in the extreme luminosity bins), P* is lognor- 
mally distributed with a spread of a = 0.4 dex. The mean 
value of P* for these galaxies is ~ —0.5, indicating that 
at the current rate of SF, it would take 3 Hubble times to 
form a quantity of stars sufficient to produce the current 
luminosity (equivalent to a Scalo 'b' parameter of ^ 0.3). 



3.2.1 Gas and SSFR in the Local Volume 

In order to be able to directly compare distributions in SFR 
and gas behaviour, the sample of galaxies must be the same 
in each case. Figure |6][ a) shows _R(Hl) plotted against Mb 
for galaxies in the Local Volume. 420 galaxies had a full 
complement of data available for analysis. Figure |6jb) plots 
SSFR against Mb in the same way. RC3 morphological types 
are colour-coded as indicated in the accompanying colour 
bar. The dashed vertical lines at Mb = -15 and Mb = -19 
indicate the boundari es separat i ng the different behavioural 
regions, identified bv lLee et al.l ()2007l ). The galaxies in this 
analysis with only upper limits on their Ha fluxes appear as 
downward- facing arrows. 

The distribution in 7?(Hl) shown inlBJa) shows the high 
luminosity outliers at low values of Mhi/M,, mirroring the 
same behaviour in SSFR. It does, however, only shows a very 
modest indication of broadening in the region Mb > —15. 

The char acteristic transit ions identified in the {EW, 
Mb} plane bv'Le e et all (j2007l ) appear in the {SSFR, Mb} 
distribution in[6jb). In particular, the distribution in SSFR 
for the least luminous (Mb > —15) galaxies is significantly 
broader than for intermediate, ( — 15 > Mb > —19) types. 
At Mb < —19, the distribution 'turns over', and spreads to 
include lower SSFR values. 



3.2.2 Gas content and SSFR of the entire sample 



The plots shown in ^3.2.11 above have the advantage of be- 
ing a very well characterised sample, with a high level of 
completeness for the volume studied. However, this benefit 
has associated weaknesses. The total volume probed by an 
11 Mpc survey is very small; approximately 5000 Mpc"^. As 
such, the sample will suffer from a paucity of rarer objects, 
such as bright, massive galaxies, and will not be a represen- 
tative sample of star formation in the 2 --^ universe. Ha 
surveys which probe greater volumes are required in order 
to provide better coverage of the upper end of the luminos- 
ity function. To this end, the larger samples described in 
i]2.1l have been included. As discussed above, surveys specif- 
ically targeting extreme or unusual systems were excluded, 
in order to provide - as much as possible - a statistically 
representative sample of star forming galaxies at 2; = 0. 

Plot El^a) plots the H I gas mass ratio against absolute 
B-band magnitude for this larger composite sample. Mb , as 
opposed to the stellar mass, has been used for the abscissa, 
in order to facil itate easier comparison with the results of 
iLee et all (|2007h . As in the llMpc sample, the distribution 
exhibits two characteristic transitions. The 'red sequence' 
transition is apparent, with the galaxies in the most lumi- 
nous bin (Mb < —19) turning off towards very low gas frac- 
tions of < 1%, and becoming dominated by early types. 
This transition luminosity corresponds to the stellar mass 
transition (3 x 10^° Mq) discussed in §1 above. At the less- 
luminous end, however, the distribution again deviates from 
the equivalent SSFR plot. Late type galaxies dominate at 
Mb > —17, having large gas fractions, typically > 50%. The 
distribution shows no sign of the broadening towards low 
values displayed by the SSFR plot, as quantified below. 

Figure[7jb) shows the distribution of SSFR with Mb for 
this total compilation of Ha data. The distributions shown 
in Figure [Tja) and Figure [7Jb) are consistent with those 
shown in Figure [6j a) and[6jb). As a qualitative check, it is 
also possible to make the above plots with rotation al velocity 
instead of Mb as the x-axis mass tracer (as in iLee et al.l 
(2007)). These are not shown for the sake of brevity, but the 
global relationships do not differ significantly from the plots 
shown. 

Plot El^b) again demonstrate s the three beh avioural 
transitions identified previously bv lLee et al] |20o3). Above 
Mb = —19, where the population is dominated by early- 
type spirals, the distribution of galaxies both broadens and 
turns over to lower SSFR values. In the intermediate lu- 
minosity bin, —15 > Mb > —19, the population is com- 
prised primarily of later type spirals and dwarfs. The mean 
SSFR is close to constant in this luminosity bin, at log 
(SSFR) ~ —10.5 yr"'^. In the lowest luminosity bin, the 
population is comprised almost entirely of late-type irregu- 
lars and dwarfs. The distribution of SSFRs broadens signif- 
icantly, occupying a range of over three orders of magnitude 
(-8.9 > log(SSFR) > -12.4). 

Figure |8] shows, for reference, the changing distribution 
of SSFR in the stellar mass plane (as opposed to using Mb as 
the abscissa). The distribution with stellar mass is very sim- 
ilar to that in the SSFR-Mb plane, including both the upper 
and lower transitions, and the interim behaviour. The 'up- 
per' transition is seen to manifest very close to t he transition 
mass, ~ 3 X 10^'^Mq previously identified (e.g. ISalim et al.l 
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Figure 6. Plot (a) - the relationship between _R{Hl) and the B-band luminosity for the Local Volume galaxies described in §1 above. 
Plot (b) - the relationship between the Specific Star Formation Rate and the B-band luminosity for the same sample. Dashed vertical 
lines at Mb = -15 and Mb = -19 indicate the boundaries separating the three SSFR behavioural regions identified by Lee et al. (2007). 
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Figure 7. Plot (a) - the relationship between the iJ(Hl) and the B-band luminosity for the complete Ha sample described in §1 above. 
Plot (b) - the relationship between the Specific Star Formation Rate and the B-band luminosity for the same sample. Dashed vertical 
lines at Mb = -15 and Mb = -19 indicate the boundaries separating the three SSFR behavioural regions identified by Lee et al. (2007). 
Least-squares best fitting power laws fitted to the distributions in (a) and (b) have indices (—0.59 ±0.02) and (—0.25 it 0.03) respectively. 



I2OO7I : ISchiminovich et "al]|2007h . and for reference a vertical 
line has been drawn at this m ass. The line definin R the ' star 
forming sequence', defined bv lSchiminovich et al.l (|2007l ) as 

log SFR/M, = -0.36 log M* - 6.4 

has been overplotted. Figure [8] is consis tent with the hig h 
mass turnover results found recently (e.g. lSalim et al.ll2007^ ■ 
which are generally found using larger samples biased to- 
wards more massive galaxies. 

The changing distributions of SSFR and Hi gas mass 
ratios are shown in more detail in Figures [9land [T0l The data 
are divided into four luminosity bins, taking the upper and 
lower marked regimes as two separate bins, and subdividing 
the middle ( — 15 > Mb > — 19) bin into two equal sections. 
In Figure m the distributions of H I gas mass ratio in the four 
bins are shown as histograms. In Figure [9l the equivalent 
histogram is shown for the SSFR distribution. The legends 



display the standard deviation of the data in the respective 
bins. 

Figure |9] displays the high-luminosity broadening in 
-R(IIl), and the lower values of a in the intermediate bins, 
mirroring similar behaviours in the respective bins of Fig- 
ure llOl There is significant difference, however, in the lowest 
luminosity bin. Whereas the distribution of SSFRs in this 
region becomes much broader than any other region, the 
distribution of Hi gas mass ratios is tighter in this region 
than elsewhere. The distributions in Hi gas mass ratio are 
also negatively skewed, as opposed to the more symmetrical 
distributions in SSFR. 

Figure [10] shows the SSFR distribution in the most lu- 
minous class to be broader than the two intermediate classes, 
corresponding to the already-noted 'turnover' towards low 
SSFRs. The intermediate bins are lognormally distributed, 
with lower values of a than either extreme bin. The lowest 
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Figure 8. The relationship between the Specific Star Formation 
Rate and stellar mass for the full compilation. The dashed vertical 
line indicates the transition mass (identified previously at ~ 3 X 
10^" Mq). The dashed sloped line represents the 'SF Sequence' 
defined by Schiminovich et al. (2007) 



luminosity bin, containing the dwarf galaxies with anoma- 
lous SF behaviour, has the highest of all standard deviations, 
and a highly non-lognormal distribution. 

Another significant difference between Figures [TOl and [9l 
is the evolution of the mean. Both the mean SSFR and the 
mean Hi mass ratio decrease with luminosity, but by sig- 
nificantly different amounts. The difference in mean SSFR 
between the two most extreme bins is 0.72 dex; this is con- 
trasted with a 1.37 dex difference in mean i?(Hl) over the 
same range. This can be more explicitly seen in the dif- 
ference between slopes for the best fitting power law: the 
i?(H l) distribution in Figure[7ja) has a power law fit of index 
A'' = -0.59 ± 0.02 (essentially identical to that in Figure O, 
whereas the distribution in SSFR is significantly shallower, 
at iV = -0.25 ± 0.03. It is worth noting in Figures l9l and [TOl 
that for all bins save the least luminous, the value of o-(Hl) 
is systematically greater than (j(SSFR). This likely results 
from the slopes of the distributions artificially introducing 
scatter across the relatively wide (2 magnitudes, 0.8 dex in 
luminosity) bins. This effect can be removed by subtracting 
0.235 from o-(E[l) and 0.1 from a(SSFR), which removes the 
mean effect of the slopes discussed below. When this is taken 
into account, cr(Hl) is lower than a(SSFR) in all bins. 

This variation in slope can be seen as a change in the 
mean Hi consumption timescale (as discussed by, f or ex- 
ample, [Robert s|[l969l : ISandage|[l98g : iKennicuttllTgSSI ): with 
increasing luminosity, the H I content drops off faster than 
the SFR, resulting in shorter H I consumption timescales for 
more luminous galaxies. This is shown explicitly in Figure 
111! where the mean Hi consumption time decreases mono- 
tonically with increasing luminosity (upward facing arrows 
denote galaxies with only upper limits on their SFRs). An 
interesting feature of Figure [TT] is the lower limit - with the 
exception of one galaxy (NGC 2681), all points have Hi con- 
sumption times >100 Myr, which is (approximately) one dy- 
namical time for a typical galaxy. This lower cut-off effect is 
caused by the SFR ceiling, found to be roughly equal to the 
gas mass divided by the minimum time to assemble the gas 
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Figure 9. Histogram showing the distribution of iJ(H l) for the 
complete Ha sample shown in Figure[7la) above. The sample has 
been divided into four luminosity bins as shown in the diagram. 
The standard deviation (cr) and the mean (<x>) are also shown; 
the position of the mean is marked on the plot with a black line. 
The value of a has been corrected to account for the slope of the 
distribution as detailed in the text. 



- i.e. the free-fall collapse time. Naturally, therefore, there 
should be a rough lower limit to the consumption time equal 
to the dynamical time. 

The H I consumption timescales present in some low lu- 
minosity dwarf galaxies are in some cases extremely high - 
for example, the dwarf galaxy UGC 07298 (Mb = -12.27) 
has a Hi consumption time of 2.7 x lO'^^ years. These long 
consumption times likely result from a combination of a very 
low SFR and an extended H I disc, where a large fraction of 
the Hi disc, though recorded by the single dish measure- 
ments, is not actively available for star formation. These 
two factors can conspire to produce an extremely long con- 
sumption time. 



3.2.3 Correlation of SFR with different ISM components 

It is a matter of ongoing debate whethe r the S F R cor - 
relates better with Hi or H2 density: iBosellil l|l994 ). 
iDeharveng et al.l (|l994l ). and lKennicufr3 (|l998bl ) find a sig- 
nificantly better correlation with the atomic gas density. 
lKennicut3 (|l998bl ) suggest that the large proportion of the 
ISM composed of H I is responsible for the better correla- 
tion with Hi while variations in the CO / H2 conversion 
are respo nsible for the poor co rrelation with molecular gas. 
However, IWong fc Bli"t3 (|2002l ) investigated the correlation 
between SFR density and gas density and found that Esfr 
correlated much better with than with Ehi, to the ex- 
tent that their H I data were inconsistent with a Schmidt law, 
while their CO data showed excellent correlation. (They sug- 
gest that this is due to Hi self shielding and converting to 
H2 at higher column densities, as described by, for example, 
iFederman eraLlll979l : IShava fc Federman|[l987l .) Q 



^ It is often the case that global studios find a superior connection 
with H I while resolved studies of individual galaxies find H2 to 
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Figure 10. Histogram showing the distribution of SSFR for the 
complete Ha sample shown in FigureOb) above. The sample has 
been divided into four luminosity bins as shown in the diagram. 
The standard deviation (cr) and the mean (<x>) are also shown; 
the position of the mean is marked on the plot with a black line. 
The value of a has been corrected to account for the slope of the 
distribution as detailed in the text. 




Figure 12. The variation of the Pearson product-moment cor- 
relation coefficient with absolute B-band magnitude, for correla- 
tions between different gas mass ratios (H I, H2 traced via CO, and 
total gas) and SSFR. Error bars were estimated using a Monte 
Carlo bootstrapping method. Due to the low numbers of low- 
luminosity galaxies, galaxies with Mb > — 13 were considered as 
one group, as were galaxies with Mb < —21. 
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Figure 11. Characteristic Hi consumption timescale (defined as 
M(Hl) / SFR and expressed in Gyr). The dashed horizontal line 
represents the Hubble time. These timescales were derived using 
a Salpeter I MF - timescales sho uld be increased by 1.44 for a 
Kroupa IMF llKroupa et al.lll993h 

In order to more rigourously quantify the relation be- 
tween gas content and SSFR, a correlation function will 
be used: being the standard test of correlation, the Pear- 
son product-moment correlation coefficient (PPMCC) will 
be usec|f|. This is defined as 

correlate better: see iFumagalli et al.l | |2009|) for a recent discussion 
of this interesting problem. 

^ Is is also possible to use a non-parametric statistic, such as 
Spearman's Rank - both measures of correlation give virtually 
identical results. 



E - (E ^0 E yf - (E y^r 




where r^y is the correlation coefficient, a is the standard 
deviation, and n is the number of points. The three di- 
mensional parameter space (= R{gas), SSFR, Mb) has been 
sliced into 1 mag bins along the Mb 'axis', and the resulting 
correlation between R{gas) and SSFR calculated for mem- 
bers of each bin separately. This way, we explicitly calculate 
the variation of the dependence of SSFR on i?(gas), as Mb 
changes. 

Due to the low numbers available at the extreme ends, 
galaxies with Mb > —13 were considered to be a single 
group, as were galaxies with Mb < —21. Figure \n\ shows 
the variation of the PPMCC (between SSFR and gas mass 
ratio) with absolute magnitude, for the three definitions of 
gas ratio discussed above. The red line shows the correlation 
between 7?H I and SSFR - in order to include as much data 
as possible, the sample drawn upon for this calculation was 
the full sample, shown in Figure [T] The green and blue lines 
show the value for 7?(H2), and 7?(total) respectively. Due to 
the lack of CO data at the lower luminosities (see Figure Q, 
the plot above has been truncated below Mb = — 16 for the 
molecular and total gas data, la errors - or more precisely 
robustness measures, being as they were estimated using a 
Monte Carlo bootstrapping method - were typically < ±0.1, 
and are shown for the H I data. For datasets smaller than ~ 
20 points, the bootstrap method of error estimation becomes 
unreliable: for this reason, bootstrap error bars were not 
calculated for molecular and total correlations. 
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4 DISCUSSION 
4.1 Error Budget 

Each observational parameter will come with some associ- 
ated error, as will each variable conversion. As a result, it is 
imperative to assess how much scatter in the relationships 
will be driven by these errors, and how much results from 
innate parametrical variance. 

Photometric errors are, in general, small, being typi- 
cally ± 0.1 mag (with the exception of the faintest dwarf 
galaxies, which can have photometric errors of ^^0.3 mag), 
as are errors on Ha and Hi fluxes. Ha extinction correc- 
tion errors are also, in general, small, though the scatter 
in the statistical relationship does increase as a function of 
luminosity. While the scatter for the least luminous dwarf 
galaxies (which have very small levels of extinction) is neg- 
ligable (~8%), the scatter at the luminous end (Mb < —19) 
is ^^40%. The scatter between these extremes is ~20%. 

By far the dominant source of error in this analysis is the 
adoption of an assumed colour in order to calculate the stel- 
lar mass. There is a a of ~ 0.3 dex in the relationship shown 
in Figure [2] which far outweighs the other error sources. As 
noted previously, however, the demographic results of the 
analysis are identical when only galaxies with true colours 
are considered; the colour adoption method will generally 
cause scatter to reduce, and indeed the most extreme points 
in all the plots in this analysis are derived from galaxies with 
true colours. We can say, therefore, that the predominant 
cause of the scatter in the diagrams is true cosmic scatter, 
and an estimate of the error in gas and SSFR parameters 
would be ±0.3 dex. 



4.2 Gas Distribution 

Figure |3] shows the distribution of the normalised H I con- 
tent for a very large sample of low redshift galaxies, across 
the entire l uminosity function. It rep roduces the long known 
result (e.g. [Roberts fc Havneslll994 l that earlier type galax- 
ies, at the more luminous end of the spectrum, have far 
less gas per unit stellar mass than their less massive com- 
panions, having long ago converted the ISM material into 
stars. The least massive galaxies, on the other hand, have 
yet to exhaust - or indeed fully tap - the H I available in the 
ISM. This is, of course, in accordance with the downsizing 
paradigm, whereby massive galaxies finish their star form- 
ing era at higher redshifts, whereas the mean star formation 
rate of dwarf galaxies has changed little over cosmic time. 

This investigation is somewhat hampered by the 
paucity of massive very early type galaxies, which comprise 
the classic red sequence. It is possible to discern, however, 
a small population of massive galaxies (Mb > —19), which 
'turn over' towards very low H I mass ratios of < 1%. This 
suggests that the red-blue sequence divide is at least con- 
nected to the availability of H I in the ISM. 

The three plots in Figure|4]show the variation in i?(gas) 
with the species of gas being considered. Panel (a) confirms 
that the H I properties of this restricted sample (only galax- 
ies with both Hi and H2 data are included here) are con- 
sistent with the larger sample shown in Figure [3] Panel (b) 
shows the equivalent plot for the normalised molecular gas 
content. Immediately obvious is the fact that there is no 



strong trend in -R(H2) with Mb, as there is in H I; the dwarf 
galaxies have normalised gas contents similarly low as the 
more massive spirals. This is, of course, an affirmation of 
the result that dwarf galaxies are Hi dominated, with the 
molecular contribution to their ISM being small, and has 
been noted previously (e.g. lYoung fc Scovillelll99"lh . It does, 
however, rely on the assumption of the infallibility of the 
CO-H2 ratio, which some authors have theorised to become 
unreliable in low mass (hence low metallicity) systems. Work 
on calibrating the co nversion factor to low metallic ity sys- 
tems is ongoing (e.g. lObreschkow fc Rawlings|[2009h . Even 
with a more sophisticated calibration, however, it is unlikely 
that the molecular component of the lowest mass galaxies is 
comparable to the Hi content. 

More interesting is the existence of an apparent max- 
imum H2 content of approximately 30-50% (log 7?(H2) ~ 
0), where no such maximum is seen in the Hi distribu- 
tion. Hi and H2 display very different behavioural features 
when spatially mapped - while the H2 distribution tends 
to be centrally peaked, falling off with radius, Hi typi- 
cally has a much more flat profile, often having a con- 
stant surf ace density out to many times the optical ra- 
dius (e.g. 'Morris fc Rickard' Il982l: iKenney fc Youngj Il989l : 
lYoung fc S covillc 1991; Beg um et al.ll2008l ). Single dish in- 
tegrated measurements of the Hi then, measure the total 
neutral hydrogen content of a galaxy, not all of which may 
be actively involved in the star formation process. It is pos- 
sible that many galaxies (particularly dwarf galaxies which 
lack the ability to gravitationally collapse them) have large 
Hi reserves inhabiting their outer discs, which is unavail- 
able for star formation: this would lead to the extremely 
high Hi contents (log _R(Hl) > 1) seen in Figures [3] and 21 
By contrast, global integrated measurements of the molec- 
ular gas mass only measure the dense circumnuclear gas, 
which suffers from this apparent maximum. This limit may 
result from a gravitational threshold effect, whereby hav- 
ing a large reserve of H2 in the ISM without concomitant 
star formation is simply not a stable situation for a galaxy. 
Alternatively, it may be that photodissociation effects are 
the cause of the limit, whereby the interstellar UV radiation 
held pho todissociates any large agglomeratio ns of molecular 
gas (e.g. ISchavell2004l : iKrumholz et al.|[2009l ). The origin of 
this threshold is not fully understood, however, and a full 
exposition of the various theories certainly lies beyond the 
scope of this work. 

Panel (c) shows the distribution of the normalised total 
gas content, calculated by adding the Hi and H2 compo- 
nents and applying a standard helium correction. Unsur- 
prisingly, at low masses, this plot is essentially identical to 
the equivalent H I plot, as the galaxies in question are highly 
H I dominated. As the galaxies increase in luminosity, the H I 
component decreases while the H2 gas mass ratio shows no 
significant trend with Mb, with the result that galaxies be- 
come more H2 dominated with increasing mass. The distri- 
bution also tightens compared to the Hi locus in panel (a). 
In addition, the gas-poor outliers are much less apparent, 
although this may result from the uncertainty introduced 
by the CO upper limits - in the absence of conclusive CO 
masses for these galaxies, it cannot be said with certainty 
that the downturn is absent. 



ISM & Star Formation Demographics 13 



4.3 Star Formation Properties and Modality 

As explained in t )3.2l above, Plot (b) of both Figures [6] 
and [7| displays e ssentially the same trends as identified by 
iLee et alj |20o3), particularly with respect to the two tran- 
sitions, at Mb < —19 and Mb > —15. These transitions 
are also apparent when considering the absolute SFR (Fig- 
ure [5} , with many galaxies in these two magnitude classes 
falling outside of a 'main sequence' of star forming galaxies. 
We address these two transitions in turn. 

At the upper end of the luminosity function, the Hq 
based SSFRs exhibit the characteristic turnover towards low 
values (< 10~^^ yr"^), indicating that their star formation 
is being suppressed. It is, however, important to note that 
this turn over does not manifest in the distributions of the 
total (atomic plus molecular) gas ratio - which as Figure [121 
shows is better correlated with SSFR at the high mass end 
than H I. It is likely that this paradoxical situation is caused 
by the fact that the number of quenched galaxies is small 
(~ 25) compared to the number of galaxies in the bulk of the 
distribution. Figure [12] shows that for the majority of high 
luminosity galaxies (making up the bulk of the distribution), 
their SFR is being driven by the molecular component of 
the ISM. However, the lack of a break in the molecular (and 
total) gas distribution (Figure [H panels b and c) suggests 
that while the main sequence of star formation is driven by 
the molecular component, quenching events that suppress 
star formation primarily affect the Hi content, which leads 
to the mirrored 'turnovers', seen in both SSFR and fni- 



4-3.1 Star Formation in Dwarf Galaxies 

We move now to the likely causes of the large spread in the 
specific star formation rate distribution of the low luminos- 
ity dwarf galaxies. The possible causes of this behaviour can 
be categorised under two main headings; firstly, the vari- 
ous physical mechanisms that could cause the star forma- 
tion properties of dwarf galaxies to not only exhibit a large 
spread, but also be decoupled from their gas contents. The 
second possibility is that at low luminosities, the star forma- 
tion rate is not being faithfully traced by the Ha luminosity, 
leading to a mischaracterisation of the SFR. These will be 
addressed in turn. 

There has been much work dedicate d to the processes 
involved in star formation at low masses (jSearle et alJ 
iGerola et al? '1980*; 'Pekel fcSil^ 1 19861 : iHunter et alT 
and Mac Low & Fcrrara 1 9991 . to name just a few). Much 
of this work has centred on the role of kinematical dis- 
turbances: being as they are small systems, dwarf galax- 
ies are inherently more susceptible than their more mas- 
sive counterparts to various types of disruption. Their shal- 
lower potential wells leave them open to disturbances result- 
ing from both internal and external events and processes 
which could act t o change the star formation behaviour. 
ISearle et al.l (|l973l ) first noticed that dwarf galaxies exhib- 
ited a distinctly different mode of star formation compared 
to ' normal' spirals, bein g dominated by a 'bursting' mode, 
and lGerola et al.l l|l980l) argued that this bursting mode re- 
sults from internal feedback processes; supernovae (the rate 
of which is proportional to SFR) in low mass galaxies have 
a large chance of disrupting the ISM so that SF is quenched. 
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until it can re-collapse enough to trigger another star form- 
ing event. 

IStinson et ah ran closed box smooth particle 

hydrodynamical simulations of low mass galaxies (3.18 x 
10* Mq - 8.6 X 10^ M0), finding that the SFR in these low 
mass galaxies exhibited periodic fluctuations (which they 
termed 'breathing'), and that the amplitude of these fluctu- 
ations varied inversely with the size of the galaxy. However, 
it is unclear whether these processes are sufficient to explain 
the range of behaviours displa yed by the dwarf ga laxies in 
our sample. The simulations of lStinson et al.l (|2007l ). for ex- 
ample, show only a range of ~ 1 dex in SFR (ranging from 
1 X 10"'' Mq yr"^ to 2 x 10"^ Mq yr"^) in the lowest mass 
galaxy considered - insufficient to explain the full 3 dex 
range in SFR (and SSFR) in Figure [5] and [3 

The possibility must then be considered that the true 
SFR in dwarf galaxies is not being reliably measured by 
standard linear conversions of the Ha luminosity. A number 
of possible causes have been discussed in previous work, with 
recent attention being focussed on stochasticity in the for- 
mation of high mass, ionizing stars, and the potential of sys- 
tematic variations in t he IMF (e.g., Pfiamm-Altcnbur g et al.l 
l2007l : lLee etHI (|2009l )). In terms of stochasticity, the upper 
end of the IMF will only be sparsely sampled in galaxies 
with low SFRs. If the mean time between formation of an 
O star is greater than the mean lifetime of an O star (~5 
Myr), then there will be times where the galaxy does not 
host an O star, but lower mass stars are being formed. This 
will result in episodic Ha emission, unrelated to the true un- 
derlying SFR, or other physical properties such as the gas or 
stellar mass. Beyond pure statistical effects, there has also 
been recent debate on whether the IMF in dwarf and low 
surface brightness galaxies could be systematically deficient 
in high mass stars. 

It may be that some combination of such effects are re- 
quired in order to explain our results in the low mass regime. 
The low mass group contains many galaxies with high SS- 
FRs, which are rapidly forming gas into stars, along with the 
galaxies with apparently quenched SSFRs. In Figure 5, for 
example, the low mass group contains both galaxies on the 
'main sequence' and galaxies with extremely low SFRs of < 
lO~^M0yr-^ . In one possible scenario, the effect of vari- 
ations in the star formation histories of low mass galaxies 
on the Ha luminosity are amplified by stochasticity in the 
formation of high mass stars. During active periods of star 
formation, the Ha luminosity is more likely to accurately 
trace the SFR, whereas during relatively quiescent periods, 
when the SFR drops below a certain threshold (~10~^ Mq 
yr"'^), the Ha emission drops dramatically, creating the illu- 
sion of more quenching than actually exists. Such issues are 
discussed in greater detail by Lee et al. (2009) and Tremonti 
et al. (2009, in preparation). 



4.3.2 The ISM-SFR connection 

Figure [T2l as explained in the text, shows the variation of the 
PPMCC with absolute magnitude for the three definitions 
of 7? (gas). Hi data being available for many more galaxies, 
the correlation between normalised Hi and SFR has been 
shown for the entire range of luminosities - the correlation 
between normalised H2 and SFR, and total gas content and 
SFR has been shown where data are available. 
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The correlation between H I and SFR is essentially the 
same for all bins with Mb < —15 - it then abruptly drops at 
the lowest magnitudes - indicative of the extreme broad- 
ening in SSFR (discussed above), which is not apparent 
in the Hi distribution. In high mass galaxies Mb < —19, 
the correlation with the molecular component is found to 
be better than the neutral component. For galaxies with 
Mb < —19, however, the correlation with the molecular 
component drops to negligible values, with the correlation 
with the ij (total) being determined entirely by (and there- 
fore approximately equal to) the value for H I. The best cor- 
relation for the most luminous bin of galaxies (Mb < —19) 
is found to be with the total gas ratio - the sum of H I and 
H2, corrected for He. 



5 CONCLUSIONS 

In this work, we have presented an analysis of the gas con- 
tent (both molecular and neutral) and Hq based star forma- 
tion rate of a sample of galaxies representative of the over- 
all galaxy population in the local universe. Several archival 
datasets and surveys were drawn upon in order to compile 
our sample, including the Ha survey of the Local Volume 
l|Kennicutt et al.1l2008l ) which provides statistically complete 
Ha coverage of field galaxies to masses below that of the 
SMC for the first time, allowing a full characterisation of the 
star formation behaviour in all types of local galaxies. The 
final sample of around a thousand local galaxies is treated in 
a homogenous and consistant manner, allowing star forma- 
tion rates, gas contents, and stellar masses to be compared 
in a uniform manner over a mass range of over 5 orders of 
magnitude. Our main conclusions are as follows: 

• Hi gas mass to stellar mass ratios (_R(Hl), defined as 
Mhi/M«) were analysed for a large sample of 2 ~ galaxies. 
7? (Hi) was found to decrease with increasing luminosity, in 
line with previous results. A small population of galaxies at 
the high mass end have anomalously low 7?(Hl), resulting 
in the distribution having a population of distinct outliers 
at Mb ~ -19. 

• The molecular gas mass ratio (-R(H2) = Mh2/M*) 
does not display the monotonic decline with increasing 
luminosity exhibited by the neutral gas. Moreover, there 
appears to a a maximal value of the molecular gas to stellar 
mass ratio (which does not exist in Hi) at ~30-50% of the 
stellar mass. The total gas distribution shows a decline with 
increasing luminosity, but no outliers; those galaxies with 
the lowest Hi masses have the largest H2 corrections. 

• The rate of star formation traces a monotonically 
increasing sequence with luminosity, with some distinct 
exceptions. At the low mass end, many galaxies appear 
to have extremally low SFRs (if the standard linear Ha 
luminosity - SFR relationship is to be accepted), falling 
off a 'main sequence' of galaxies defined byl>P*>— 1. 
Intermediate luminosity galaxies fall within this sequence, 
almost without exception. At the high luminosity end, in 
addition to a slow and general trend towards lower SFRs 
(and P* values), there is a distinct population of galaxies 
with 'quenched' SF, with extremely low SFRs. 



• When directly comparing gas mass ratios with SSFRs 
(so both SFR and Mhi normalised to the stellar mass) for 
identical samples of galaxies, it is clear that the population 
of low mass galaxies with a wide spread in SSFRs is not 
matched by an equivalent shift to a larger range of 7?(Hl) 
- the dwarf galaxies with extremely low SFRs do not, 
apparently, suffer from a dearth of fuel for SF in their 
ISM. At the high mass end, the SSFR matches the Hi 
distribution well, suggesting that the lowering of the SFR in 
massive galaxies is intimately connected to the Hi content. 

• Although our restricted samples preclude a full analysis 
of the origins of the red / blue sequence bimodality, it must 
be remarked that the results presented herein are consistent 
with the finding that the transition between the red and 
blue sequences is, at the very least, intimately connected 
with changes in the available gas supply, and may in fact 
be primarily driven by the amount of gas available to fuel 
star formation. 
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